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Abstract. Quantum metrology offers an enhanced performance in experiments such 
as gravitational wave-detection, magnetometry or atomic clocks frequency calibration. 
The enhancement, however, requires a delicate tuning of relevant quantum features 
such as entanglement or squeezing. For any practical application the inevitable 
impact of decoherence needs to be taken into account in order to correctly quantify 
the ultimate attainable gain in precision. We compare the applicability and the 
effectiveness of various methods of calculating the ultimate precision bounds resulting 
from the presence of decoherence. This allows us to put a number of seemingly 
unrelated concepts into a common framework and arrive at an explicit hierarchy of 
quantum metrological methods in terms of the tightness of the bounds they provide. In 
particular, we show a way to extend the techniques originally proposed in Demkowicz- 
Dobrzanski R, Kolodynski J and Gu^a M 2012 Nat. Commun. 3 1063, so that they 
can be efficiently applied not only in the asymptotic but also in the finite-number of 
particles regime. As a result we obtain a simple and direct method yielding bounds 
that interpolate between the quantum enhanced scaling characteristic for small number 
of particles and the asymptotic regime, where quantum enhancement amounts to a 
constant factor improvement. Methods are applied to numerous models including 
noisy phase and frequency estimation as well as the estimation of the decoherence 
strength itself. 



PACS numbers: 03.67.-a, 03.65.Yz, 03.65.Ud, 03.65.Ta, 42.50.St 
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1. Introduction 



Quantum enhanced metrology has recently enjoyed a great success at experimental level 
leading to new results in atomic spectroscopy magnetometry p-fTf and optical 

interferometry with prominent achievements in gravitational 



waves sensing 



12 



As predicted in the incipient theoretical results of |!3l-ll7|. a physical parameter unitarily 
encoded into a quantum system — probe — consisting of N entangled, non- inter acting 
particles (atoms or photons) can be extracted with precision that is limited by the 
quantum-mechanical uncertainty relations and not the more restrictive central limit 
theorem of classical statistics. Hence, the uncertainty in reconstructing the encoded 
parameter, such as e.g. optical phase delay or frequency difference, can in principle be 
proportional to 1/N, the so-called Heisenberg Limit (HL), rather than l/y/N, commonly 
referred to as the Standard Quantum Limit (SQL) or the shot (projection) noise. 
However, this dramatic scaling improvement can be illusive, as both the experimental 
results and theoretical toy-models have indicated that achieving the ideal HL is a hard 
task owing to the strong destructive impact of imperfections, which should be always 
accounted for in realistic scenarios. 

An important question that has been considerably addressed by many researchers 



18H32fl reads: how and to what extent can the noise effects be compensated in 



quantum metrological setups? In the case of atomic spectroscopy, it has already been 
indicated in 24|, |25[ that the effects of uncorrelated noise, independently affecting 
the atoms within a probe, have a dramatic impact on quantum protocols - - most 
likely restricting the ultimate precision scaling to become SQL-like for high enough 
N, so that the quantum enhancement is asymptotically limited to a multiplicative 
constant factor. In optical interferometry, photonic loss is the main obstacle to practical 
implementations of quantum enhanced protocols 27H30] and the asymptotic SQL-like 
scaling is again inevitable, as proven in 
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|. Similarly to the atomic case, the 
asymptotic improvement constant becomes an essential feature which determines the 
achievable precision for high N. Following the above exemplary models, methods of 
quantifying the asymptotic quantum enhancement have been proposed for arbitrary 
kinds of probes with decoherence present 33M39|. Recently, general procedures have 
been established that are capable of deriving practical bounds on ultimate achievable 
precision in realistic quantum metrological setups 38|, |39|. In particular, in the case of 
uncorrelated noise it has been demonstrated that the sole analysis of the evolution 
of a single particle often leads to surprisingly informative bounds on the precision 



achievable with arbitrarily entangled multi-particle inputs 39] - see Figure [T] for a 
summary of the single particle (single channel) methods investigated further on in this 
paper. For completeness, it should be noted that there are some specific metrological 
models with noise, in which asymptotic scaling power enhancement is nevertheless 
possible 22|, |23|, |40| . Still, applicability of these models is limited, since in any practical 
implementation the noise types considered will always be accompanied by some generic 
decoherence processes for which the constant factor bound on the maximal quantum 
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Figure 1. Quantum metrology and the single channel methods 

(a) General scheme for quantum-enhanced metrology with uncorrelated noise. The N 
particles within the probe in a quantum state ipF evolve and decohere independently 
while sensing an unknown parameter (p (e.g. phase). An estimate p of the parameter 
is inferred from a measurement result on the final state of the probe . 

(b) Precision bounds from single channel analysis. The ultimate precision is bounded 
by the 1/vN scaling (SQL), if for small variations bp around ip the channel can 
be expressed as a parameter-independent map <f> that is also fed a classical, diagonal 
state Pu>, classical simulation (CS), or a general, quantum state a v , quantum simulation 
(QS), varying smoothly with p. Still, for all such channels and more, the tightest bound 
on precision is obtained by the channel extension (CE) method, in which the map A v 
is replaced by its extension A v ® X. 



enhancement will force the asymptotic precision scaling to be SQL-like. 

In this paper, we provide new insight into relations between seemingly unrelated 
methods used for derivations of various quantum metrological bounds and order them 
with respect to their predictive power. Firstly, focussing on the geometric intuitive 
method of channel classical simulation (CS) introduced in 36|, |39(, we prove that its 
applicability is equivalent to the approach proposed in j3?]|. Moreover, we show that 
the criterion for classical simulability of a channel can be generalized to a quantum 
simulability (QS ) condition (36|, which coincides with the channel programmability 
postulate of [34]. Although the idea of QS allows to prove the asymptotic SQL-like 
behaviour for a wider range of decoherence models, we demonstrate that the channel 
extent (CE) method of encompasses all CS, 0, QS and approaches 
providing more stringent bounds on precision. We also comment on the problem of 
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the tightness of these bounds, as they are guaranteed to be saturable only for channels, 
for which the estimation task cannot be improved by allowing for an additional ancilla 
(extension, idler mode) entangled with the input state. Classification of such channels 
is an open problem of current research in channel estimation 4jJ-|44j and akin channel 
discrimination 45 53| theory. The graphical interpretation of the CS, QS and CE 
methods is presented in Figure [TJa. 

Most importantly, we go beyond the results of 39| and show that the CE method 
may be applied not only in the asymptotic regime but also when dealing with finite 
probes of N particles. Similarly to the asymptotic case, it corresponds then to an 
optimisation procedure over Kraus representations of a given channel that can be 
recast into an efficiently solvable semi-definite programming task. We apply our 
results to phase/frequency estimation with various noise models including: dephasing, 
depolarization, loss and spontaneous emission; restricting ourselves to the cases of noise 
commuting with the parameterized unitary part of the evolution. This assumption 
makes the analysis more transparent, but is not indispensable, since our methods may 
be effectively employed for any single particle evolution described by a g eneral Lindblad 
equation |54J reshaped into the corresponding Kraus representation (55|. What is more, 
as our finite- iV CE method applies to models for which its asymptotic version fails, it 
can be used to up per -bound the asymptotic scaling for channels surpassing the SQL. 
As an example, in 4o|, our new method has been already utilized to predict asymptotic 
super-classical scaling for a channel with a non- commutative noise. Additionally, in 
order to stress the generality of the methods, in the final section of this paper we show 
that they can be applied not only to noisy unitary parameter estimation tasks but also 
to ones in which the decoherence strength itself is estimated. Finally, we should also 
clarify that noise correlation and memory effects 56|, |57fl are beyond the scope of this 
paper and while the non-Markovian effects, provided they affect each of the particles 
independently, might be analyzed using the tools presented, correlations between the 
decoherence processes affecting the particles does not fit well into the framework here 
advocated. 



2. Quantum Fisher Information 

2.1. Classical Cramer-Rao bound 

Let us assume that after measuring a physical system an outcome is obtained that can 
be represented by a random variable X distributed with some probability distribution 
Py(X). If the system is classically described, all its properties can be simultaneously 
determined, so that X can in principle be multidimensional and contain as much 
information about the system as allowed by the available resources. The estimation 
task corresponds then to determining with highest precision the quantity tp based on the 
observed value of X. As stated by the Cramer-Rao bound (58| any unbiased strategy to 
determine the unknown parameter after repeating the procedure k times, must provide 
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an estimate <p with uncertainty that is lower bounded by 

> J— , where F^} = fdx (1) 

is the (classical) Fisher Information (FI) 0. 

The 1 / dependence in ([1]) is a consequence of the central limit theorem and 
the fact that the k procedures are independent. This manifests itself by the additivity 
property of the FI, i.e. F c i[p^j = k F^p^] for X k . Equation ([1]) shows that the FI 
is a local quantity containing information about infinitesimal variations of ip. That is 
why, FI is designed to be used in the so called local estimation approach in which small 
parameter fluctuations are to be sensed. This small deviations regime may always be 
reached after many procedure repetitions (k — > oo) and in this limit the Cramer-Rao 
bound is known to be saturable via e.g. max-likelihood estimation schemes |58| . 

2.2. Quantum Cramer-Rao bound 

In a quantum estimation scenario, the parameter ip is encoded in a quantum state p v . 
A general measurement, mathematically represented by the elements of the positive 
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60 , is 



operator valued measure (POVMs) M x that satisfy M x > 0, fdx M x = I [jj 
performed yielding outcome statistics p v {X) = Tr {p v Mx}- Establishing the optimal 
estimation strategy corresponds then not only to the correct interpretation of the 
measurement results, but also to a non-trivial optimization over the class of all POVMs 
to find the measurement scheme maximizing the precision. In this case the Quantum 
Cramer- Rao bound can be derived 61 63|, which is independent of the choice of the 
POVMs and solely determined by the dependence of the output state on the estimated 
parameter: 

A^>^= with F q [p v ] = Tt\p v L s [ Pip ] 2 } (2) 

being now the Quantum Fisher Information (QFI). The Hermitian operator L s [p v ] is the 
so called Symmetric Logarithmic Derivative (SLD), which can be unambiguously defined 
for any state p v via the relation f> ip = \ (/^-^[pj +L s [p ip ]p^) . Then, in the eigenbasis of 
P<P = J2i^i( i P) kiMX^MI with {|e i ((/3))} i forming a complete basis (V,: < Ai < 1) 

r s r i V- 2 (eifaOl A? Ml)) i / U/ / N, f o\ 

H = L x^pJTxM ^))(^)\ , (3) 

where the sum is taken over the terms with non-vanishing denominator. QFI is additive 
on product states and in particular Fq [p® k ] = k Fq[ Pip ]. Thus, the \fk term in the 
denominator of §Z§ may be equivalently interpreted as the number of independent 
repetitions of an experiment with a state p v o r a single shot experiment with a multi- 



party state p® . Crucially, as proven in [63, |64|, there always exist a measurement 



% Throughout the paper, we depict derivatives w.r.t. the estimated parameter with an 'overdot', so 
that e.g. p v (x) = d v p v (x), p v = d v p v and K(ip) = d v K(ip). 

§ We denote by I — the identity operator and by I — the identity superoperator. 
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strategy, e.g. projection measurement in the eigenbasis of the SLD, for which bounds 
P]) and (j2J) coincide. Hence, as in the classical case, the saturability of ([2D is guaranteed, 
but again only in the k— >oo limit. 

2.3. Purification-based definition of QFI 

For pure states, p v = |?/v)(Vvl' the Q^I in ([2]) simplifies to -Fq[|^)] = 
4 y(^pifi\'ip(p^ ~ (i ) <p\'4 ) <p^\ J SI- Yet, as indicated by ([3]), for general mixed states the 

computation of QFI involves diagonalisation of p v which may be infeasible when dealing 
with large systems. That is why, it is often necessary to look for upper bounds on 
QFI that would be efficiently calculable even at the expense of saturability. For this 
purpose, definitions of QFI were proposed that do not involve computing the SLD, but 
are specified at the level of state purifications: p v = Tr E {|\l/(<p))(\l/(<p)|}. In j38| the QFI 
of any p v has been proven to be equal to the smallest QFI of its purifications \ty((p)} [jj 

F Q [p v ] = minF Q [|*(<p))] =4nun{(*( V >)|%0) " |(*MI*M}| 2 }- (4) 
Independently, in (3^) another purification-based QFI definition has been constructed: 



F Q [p v ] = 4min(^)|*H). (5) 

Despite apparent difference, Eqs. (j3|) and ([S]) are equivalent and one can prove that 
any purification minimizing one of them is likewise optimal for the other and satisfies 
the condition &(<p)^ = ^L s [p v ] (gilg \^((p)) causing the second term of (j3J) to vanish. 
Although for any suboptimal ^((p) (jlj) must provide a strictly tighter bound on QFI 
than ([S]), the latter definition, owing to its elegant form, allows for more agility in the 



minimization procedure, so that it has been efficiently utilized in 35|, |39j and is also the 
base for this paper. 

2-4- RLD-based upper bound on QFI 

On the other hand, a natural way to construct a bound on QFI and avoid the SLD 
computation is to relax the Hermiticity condition of the logarithmic derivative. If a 
non-Hermitian L[p v ] satisfying d^p^ = \ ( < p (p L[p lp \+ L[p,p]^ p^\ can be found, as proven 
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62 . |64j, an upper limit on QFI in ([2]) is obtained: -Fq[pJ < Tr |p v L[p ¥ ,]L[p v ,]t }. In 
particular, if and only if d^p^ is contained within the support of p v , one can construct the 
Right Logarithmic Derivative (RLD) by setting L[p v ] =L R [p v ] = p~ x d^p^ and formulate 
an upper bound on QFI of a simpler form: 

FqW] < -^Q LD \p<p\ = Trjp" 1 (<W}- (6) 

I We shorten the notation of functions and superoperators of pure states, so that i r [|V')] = 
andA[|^)]=A[^)(V|]. 

% See also an alternative formulation based on the convex roof formula, which is valid for unitary 
parameter estimation j(35l. [66|. 
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Although ([6]) is tight only when -^[pj = L s [p v ], it still allows to quantify precision 
well for channel estimation tasks (37j |. as described in the following section. Lastly, one 
should note that we are not considering here multi-parameter estimation schemes, for 



which the RLD may sometimes provide tighter bounds than the SLD 67J, |68 



3. Estimation of a quantum channel 

3.1. Channel QFI 

As in metrological setups the estimated parameter is encoded in the evolution of a 
system, we identify p ip = A ip [p in } as the final state of a system that started from an input 
Pi n . The preparation of pj n is controlled in order to achieve the most precise estimate 
of <p that parametrizes some general channel — a Completely Positive Trace Preserving 



(CPTP) map Ay, |59l . |60|. Although the form of A^ in general strongly depends on the 
model considered, the Quantum Cramer-Rao bound always applies, so that precision 
is upper bounded according to (J2J) with QFI ^[A^fpin]]. Furthermore, as the QFI is a 



convex quantity |69|, one should restrict oneself to pure input states when seeking the 
optimal one. Hence, as shown in Figure [2^, we define the channel QFI as the maximal 
QFI after performing the input optimisation, so that it has a concrete operational and 
application-like interpretation: 

F[K V ] =maxF Q [A,[|^)]]. (7) 

fin 

For instance, while estimating the duration of the evolution (<p = t) in an ideal, 
decoherence-free setting, the CPTP map is unitary leading to a pure channel output. 
Thus, A^fl^in)] = U t [\ip-m)} = e~ im |^in)(^in| e iHt with H being the Hamiltonian of the 
evolution. Hence, the definition (J2J) corresponds to 

T\U t \ = max F Q [e- iHt \4j in )} = 4max{(^ n |i/> in ) - ^ in \H\^ in ) 2 } = 4max A 2 H, (8) 

^in Vto V>in 

and the optimal states are the ones that maximize the Hamiltonian variance. Note 
also that in this case the Quantum Cramer-Rao bound takes the form of the time- 



energy uncertainty relation, AH -At > 1/2 70|, [71], with At being the uncertainty in the 
estimated duration. 

3.2. Purification-based definition of channel QFI 



In order to employ the definition (jSJ), we utilize the Stinespring theorem [60| and express 
the channel A^ as a unitary map U^ E on the system combined with an environment 
disregarded after the evolution. In this way, the output state and its purification 
respectively read A^l^m)] = Tr E {|^(<p))(^(v^)|} and \V(<p)) = U* E \ip ln ) ® |1), where 
|1) is an arbitrary fixed state chosen to be the first vector in the basis {|«)}[ =1 of 
the environment Hilbert space W^. By specifying the dimension r of to be equal 
to the rank of A^, we can generate all non-trivial purifications, \l/(<p), by applying 
a fictitious, possibly <p-dependent unitary rotation, vE, to the environment, so that 
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Figure 2. Channel QFI based on the output state purification 

(a) Channel QFI as the QFI of the output state maximized over all pure input states. 

(b) Channel QFI obtained from the output state purification generated by a local, 
fictitious, parameter-dependent environment rotation. 

(c) Extended channel QFI independent of the maximization over the input states. The 
environment rotation corresponds to a choice of Kraus representation of the channel. 



^(y?)^ = U^ E |^in)®|l) with U^ E = u E U^ E . Furthermore, writing the channel action in 

its Kraus representation form A^U^n)] = Y7i=i Ki(f) iV'inXV'inl Ki(ip)\ we can identify 
the Kraus operators corresponding to ^(ip) as 

(9) 



KM 



(*\uT |i> 



3=1 



where u E ^ 



U u 



\j) and Kj((p) = (j\lfz E \l) are the Kraus operators of the 
original purification ^((f). Hence, by picking an environment unitary rotation u E , 
we are, equivalently to the purification choice, specifying a Kraus representation of 
A„. Moreover, as the QFI is a local quantity, we can restrict ourselves to infinitesimal 



rotations, u 1 



-ih(ip-tp ) 



in the vicinity of the real value ip that are generated by some 
Hermitian h. Taking without loss of generality ip o = 0, we obtain the channel version of 
(jSJ) shown in Figure [2]b as 

r 

J"[AJ = 4max min (SSf(<p)\Sir(<p)) = 4max min V K ^tpf R fa)^ , (10) 

Win h Win h ^—f 

1=1 

where |i(</?)) = {Uf - ihUf)\^- m ) ® |1) and similarly K t (ip) =K i {ip) - i £J =1 h^K^ip) 
with hij = (i\ h \ j) being the elements of the generator of Kraus representation rotations 
([9]). Crucially, Figure [2j) and Equation (TIT)]) indicate that the optimal purification/ Kraus 
representation corresponds to the choice of an artificial environment that rotates locally 
with if hindering as much as possible information about the estimated parameter. 
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In order to make the reader familiar with the above formalism, we apply the 
definition (fTUj) to the previously mentioned case of the evolution duration estimation 
Wt[|^>in)]- As the evolution is unitary and r = 1, the environment /Kraus rotations 
correspond just to a phase variation uf = e~ lht with h being a real scalar. Hence, 
f fTUj) simplifies to ([8]) as expected: 

F[U t ] = 4maxmin{(V>i„|#Vin> - 2/j (^in|##in) + h 2 } = 4 max A 2 H (11) 

i/) in h i/) in 

with minimum occurring at h = (tl) in \H\ip in ) . Consistently, we would also arrive at ([S]), 
if we had used the other QFI purification-based definition (J4]) as shown in 



3.3. Extended channel QFI 

The channel QFI ([7]) is affected by any ^-variations in the form of A^, quantifying 
the distinguishability between maps A^, and A^+5^. However, any such disturbance 
may be noticeable only for input states which lead to a measurable change of the 
channel output that is at best in some "orthogonal direction". As a consequence, the 
quantity min^j. . .} in (fTUj) depends strongly on the pure input -0^, as the minimum 
occurs for Kraus operators |A'° pt ((^)}^_ 1 which fulfill the condition K op \(p) \ip- m ) = 
7jL s [Ap [|^in)]] K° p \cp) \ip- m ) required for the purification of (j3J) and (j5J) to be optimal. 
Maximization of this quantity over \ip- m ) is difficult in general, due to the impossibility 



35 



of exchanging the order of max and min in (TTO]) 

Yet, one may construct a natural upper bound on the channel QFI (jTJ) by extending 
the input space, "Kg, by an equally-large auxiliary space, Ha, which is unaffected by 
the map and measured along with the channel output (see Figured). This way, by 
employing extended input states entangled between these two spaces, l^fn*) £%s x ^A ; 
one may acquire full available information about tp imprinted by the map A^ on the 



extended output state. The analogue of ( II Op defines then the extended channel QFI |35 



'121 



J 7 ^ ® X] = 4maxmin Tr s < pf a ^ K i((p)^ K i((p) } = 4 min 

Pfn h { tt J h 

where ||. . .|| represents the operator norm. The first expression is obtained by tracing 
over the auxiliary space Ha, what leads to the maximisation over all mixed states 
pf n = TrA{|V ; fn t )(V ; fn t |}- However, this is exactly (TT0l) with the pure input state re pla ced 
by a general mixed one, in which case the order of max and min can be swapped |35| El- 
Consistently, if the optimal input state of (1T2"]) is pure, (fTUj) and (fT2"|) become equivalent 
manifesting the uselessness of entanglement between Tis and "Ha and irrelevance of the 
auxiliary space. 

Importantly, the extended channel QFI ()12p can always be efficiently evaluated 



numerically by means of semi-definite programming 39| and we show that this is a 



special case of a more general task of bounding A^-parallel channels QFI, as explained 

+ We should stress that (|12[) does not correspond to the situation of using mixed states as inputs for 
unextended channel, as mixed states never outperform pure state inputs due to convexity of the QFI. 
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in Appendix E For phase estimation schemes with relevant noise models including: 
dephasing, depolarization, loss and spontaneous emission; we determine analytically 
both (TT0|) and (I12p to verify if the use of extended, entangled input states may improve 
estimation at the single channel level. The corresponding unextended/extended channel 
QFIs are presented in Table [H whereas the optimal purifications yielding the minimum 



of f|T2|) can be found in Appendix A along with the details of the channels considered. 
The results justify that extension enhances the precision only for depolarization and 
spontaneous emission channels, for which the input states maximally entangled between 
Hs and Ha are optimal. 



3.4- RLD-based upper bound on extended channel QFI 



In 37| the applicability of the RLD-based bound ([6]) has been addressed in the context 
of channels. By defining the Choi-Jamiotkowski (C-J) representation 60| of a particular 
map Atp as = A^ ® X [|I)] with |I) = X^i=i ^ S N) ® it has been proven that the 
extended channel QFI can be further upper-bounded by 



7^®X] <J^ RLD [A^®X 



(13) 



where 



.|| is again the operator norm and is the inverse of f^A^ restricted only to 
the support of the C-J matrix. Most importantly, the bound (lT3"j) is determined solely by 
the form of A v , i.e. its C-J representation, and does not contain any extra information 
about the space of input states accepted by the map. This contrasts the definitions 
of purification-based unextended/extended channel QFIs ffTU]) / (lT2"j) and facilitates the 
analyticity of the results presented in Table [TJ On the other hand, as indicated in 
Section [531 both applicability and tightness of the RLD-based bounds are limited. The 
bound (TT3|) is valid only when f2 A is fully supported by [37| • However, we give 
an intuitive reason for this restriction by proving that this condition is equivalent (see 



Appendix B) to the notion of channel A^, being Lf -non- extremal, as introduced in 



and also revisited in the following section. This confirms that the exclusive dependence 
of (ITS]) on JIa has indeed a strong geometrical meaning. Moreover, although the RLD- 
based bounds depicted in Table [T] for the relevant channels seem to be far above the 
corresponding channel QFIs — ([7]) and ( lT2"j) . they are of great significance. The bound 
(TTB"]) is additive for any maps A^, A^ to which it applies 



37 



^[(A«®Z)®(A<, 2) ®J)] =T 



RLD 



(A^®X) 



< T 



RLD 



(A^®X 



thus 

[A£>®X]+J^[Ag>®Z] 
-NF RLD [A V ®1}. 



(14) 



Hence, it constrains not only the QFI of a single extended channel, but also restricts the 
QFI of N extended channels used in parallel to scale at most linearly with N. Crucially, 
as the extension can only improve the precision, (1141) is also a valid upper-bound on the 
QFI of N uses of an unextended channel, which asymptotic precision scaling is then 
limited to a constant factor improvement over the SQL (see Section ITT]) . 
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Table 1. Channel phase estimation sensitivity quantified via QFIs and their 
asymptotic bounds. Noise models of metrological relevance discussed in the paper 
are listed in the first column. Decoherence strength increases with a decrease of the 
r\ parameter (0 < 77 < 1) (see [Appendix A] for details). From left to right — single 
channel QFI; extended channel QFI; upper bounds on asymptotic channel QFI (fl"6|) 
in ascending order: channel extension bound (see Section |4. 1 .3|) . quantum simulation 
bound (see Section 14.1.21) , RLD-based bound (see Section I3.4|) . classical simulation 
bound (see Section [4.1.1 1) . [n.a. — not available] 



4. Estimation of N independent quantum channels 

In order to describe general metrological schemes depicted in Figure [T]a., we model the 
evolution of all particles within the probe as N identical, independent channels acting 
on a possibly entangled, pure input state of the whole probe \tp^)- The final output 
state of the probe then reads = A® iV [|^[)] yielding the N-channel QFI: 

7[AH=maxF Q [Ani<}]], (15) 

K 

which linear or quadratic dependence on N dictates respectively the SQL or HL scaling 
of precision. For example, when considering classical schemes that employ unentangled 
probes, \ip^) = <S>n=i l^m)' we are effectively dealing with N independent subsystems, 
so that J-^A®^] =jV.F[Ap] and the uncertainty of the estimate (p decreases as 1/y/N. 

4-1- SQL-like bounds on precision in the asymptotic N limit 

To investigate channels that incorporate the uncorrelated noise restricting the 
asymptotic precision scaling to SQL, we define the asymptotic channel QFI as 

^as[A,] = lim -LpL. (16) 

Thus, for such SQL-bound channels, (fT6|) is finite and -FasfAJ > J-'fAJ with equality 
indicating the optimality of classical estimation schemes. Hence, ([To]) quantifies the 
maximal quantum precision enhancement that reads 



X\K\ = lim ^ = J^4 > !■ (17) 
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However, as the computation of (116]) is generally infeasible owing to the complexity of 
QFI rising exponentially with N, one normally needs to construct an upper limit on 
the iV-channel QFI (USD, J^ bound [A®*] , from which the asymptotic form, J^ ound , may 
be constructed using (TT6"j) that upper-bounds both the A r -channel QFI and the maximal 
quantum precision enhancement: 

/ TrboundTA 1 

j-[ah < n j- a b s ound , x\k\ < si • ( 18 ) 

Methods of constructing J 7 ^ " 11 ' 1 were proposed in 39| basing on the concepts of 



channel classical simulation (CS) and channel extension (CE). As mentioned already, 
the CS method applies only to tp -non- extremal channels, for which also the RLD- 
based bound (1141) provides a valid example of J r ^° und . Yet, the notion of CS may be 



generalized to the idea of channel quantum simulation (QS) introduced in [36], in order 
to obtain an asymptotic bound applicable to a wider class of quantum maps. All these 
four approaches (J 7 ^ 3 "^ = J^ s , J-^ E , J rRLD [A v ,®X] ,J-® respectively) are presented in 
Table [T]on the right hand side of the double-line for the relevant channels. As the lossy 
and spontaneous emission interferometry cases are examples of y>extremal maps, they 
do not allow for J 7 ^ and J rRLD [A ¥ ,(8)X] to be constructed. In the case of depolarization 
channel, which is full-rank (go(| and hence ^-non-extremal, despite the lack of a simple 
geometric interpretation of its value, J rRLB [A^ ®X] proves to be tighter than J 7 ^. The 
more general QS approach not only applies also to the (93-extremal) lossy interferometry 
case, but also provides as accurate bounds as the CE method. Nevertheless, as the CE 
approach is proven to provide at least as tight bounds for the broadest class of quantum 
channels containing ones to which the other methods apply, we use it to quantify the 
maximal quantum precision enhancements (ITTj) for the channels considered (see Tabled]). 
Below, we describe alternately the CS, QS and CE methods — ordered according to 
their power and generality. 

A. 1.1. Classical Simulation method 

As introduced in H and depicted in Figured),, a channel admits a dass,cal simulation 
(CS) of itself, if for any g it is expressible in the form 

Me] = %®jv] = ZX,ii%], ( 19 ) 

i 

where = ^Y^iV^.i l e «)( e «l is a classical, diagonal density matrix in some basis and $ 
is a (/^-independent CPTP map acting on a larger input space. By defining Hi[g] = 
5>[g(g)|ej)(ej|] one obtains the second equality in (j!9p . so that it becomes evident that 
the estimated (p parametrizes only the mixing probabilities of some (^-independent 
quantum maps. Then, the A r -channel QFI ( fT5|) can be simply upper-bounded via 
jr[A® N ] < NF cl [p^i], where FJv^is the discrete version of classical FI in ([T]) and 
plays the role of J 7 as° und in dUD [36j,|39j. Moreover, as QFI is a local quantity, in order to 
construct a CS-based -7 r a J s ound valid for small deviations 8<p around a given ip, the classical 
simulation must be feasible only locally, i.e. A^g] — $^j£W ^i[o\ + 0(5ip 2 ). Therefore, 
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as proven in |36|, if the C-J representation Q,^ of a channel A v allows for parameters 
e± > such that the matrices f2n± =^a v ± e±Cl^ are positive semi-definite at a given tp, 
the channel is (p -non- extremal there and the necessary p^^ can always be found. This is 
because one can construct f2 A ^ =p !p+ Qu + + 2V,-^n_ that up to 0(5p 2 ) coincides with 
Qa v by choosing p Vi ± such that —Q^ and fi A —Cl\ v . Hence, A^g] = A v [g] + 0{8ip 2 ) 
with A^[g\ = p lfit+ Il + +p v _n_, so that F c ] |? v+ 1 = l/(e+e_) is a legal example of the 
required bound valid at ip. Furthermore, in [39|, it has been proven that for channels 
of the form A v \g] = A[Up[Q]] this two-point construction is always optimal at any ip 
when maximal possible e± are chosen Q Geometrically, imagining the convex set of all 
CPTP maps in their C-J representation that share input and output spaces of A^, the 
channels f2n ± he at the intersection points of the tangent generated by VL^ v at Vt\ and 
the boundary of the set. The J^ s ound of (TTS]), which we refer to as the asymptotic CS 
bound — J 7 ^ = l/(e™ ax e max ), is dictated then by the "distances" e™ ax of the channel 
to the boundary measured along this tangent. Although the CS approach provides 
weaker bounds than the CE method (39)], it gives an intuitive geometric explanation of 
the inevitable asymptotic SQL- like scaling of all <£>-non-extremal maps. These naturally 
include the full-rank channels 
boundary. 

4-1.2. Quantum Simulation method 



60|, which lie inside the set of CPTP maps away from its 



In [36|, a natural generalization of the channel CS has been proposed which is 
schematically presented in Figure QJ). This, so called, quantum simulation (QS) of 
a channel corresponds to expressing its action in a form similar to f)19p that reads 

A v [g] = % ® a v ] = Tt e ^{u (<? ® IW>e*eA|) (20) 

where now oy is a quantum, non-diagonal, ^-dependent density matrix, and the purified 
form on the right hand side involves both channel $ environment space E$ as well 
as o^p purification space E CT such that = TrE CT {|Vv)(Vvl}- Note that the notion of 
quantum simulability is equivalent to the channel proqrammability concept introduced 



111 



34|. Following the same argumentation as in 39j for CS, the A^-channel QFI (|T5|) 
of a quantum simulable channel can be linearly bounded as J r [A? )JV ] < A r i 7 Q[a ¥ ,], and 
therefore the asymptotic bound reads J r ^° nnd = Fq [aj . Similarly to CS, a channel may 
admit many decompositions (120]) and the optimal one must yield the lowest Fq^]. 
Therefore, without loss of generality, in the search for the optimal QS, we may take 
U in (120]) to act on the full purified system, i.e. also in E$ and E CT spaces. This 
enlarges the set of all possible QSs beyond the original ones U = U SF " t ®I E ' T and yields 
^ob ~-^q[|Vv)]i which via (jl]) cannot be smaller than Fqfa^]. In fact, (jl]) assures that 
for any QS employing a v , there exists an "enlarged" decomposition f )20|) leading to the 
same J r ^° und with being the minimal purification in (jl]). In conclusion, we may seek 
the optimal QS by analysing all possible decompositions of the form (120]) that, owing 

* Yet, it may prove optimal for channels of other type, as shown for transversal dephasing in (4ol |. 
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to the locality of the QFI, must be feasible only for small deviations 5<p from a given tp, 
so that A v [g] = Tr E$E(T {£/ (£><8>|Vv)(Vvl) ^ } + 0(5<^ 2 ). In |Appendix D[ we prove that, 



in order for the QS (120]) to be possible locally at </?, A^ of rank r must admit Kraus 
operators {Ki(ip)} r i=1 that satisfy the two conditions: 

r r 1 

i^^(^(v9) = and ^^(^(^^-FqO^)] I (21) 

t=l i=l 

Hence, by optimizing over all locally equivalent Kraus representations of A^ — the ones 
related to one another by rotations (J9]) generated by any Hermitian h — that satisfy 
constraints (|2"T|) . we may determine the asymptotic bound given by the optimal local 
QS, which we refer to as the asymptotic QS bound — J 7 ^, as follows 

J 7 ^ = min A subject to: a k = ~ I, (3^ = 0, (22) 

h 4 

where a ^ = Y2l=i K i(ip)^ K i(ip) , (3^ = i YH=i (v 9 )^ Ki (v) an d ^ nas the interpretation 
of ^as° und = F q[I^)] for a local QS of the form fl2U]). Before revisiting the CE method 
explicitly below, we should note that (122]) resembles exactly the asymptotic CE bound 



of |39| with an extra constraint forcing the operator a ^ to he proportional to identity. 
This proves that indeed the QS method can never outperform the CE approach. 

Investigating the relevant quantum maps considered in Table (TJ the QS method 
must naturally apply to dephasing and depolarization channels. These are (^-non- 
extremal (classically simulable) maps, hence their locally constructible CSs (Ti~9]) serve 
as examples of the more general QSs ( 12"U|) . Consistently, the Kraus representations 
utilized in to obtain the asymptotic CE bounds for these two channels fulfil the 
a k oc I constraint. Thus, QS is not only feasible in their case but also its asymptotic 
bound coincides with the one of the superior CE method. Significantly, also in the case 



of the lossy interferometry the optimal Kraus operators used in [39| to minimize the 
asymptotic CE bound satisfy the extra QS's constraint. This fact indicates that also 
for yj-extremal channels QS may prove to be as good as CE. However, in the case of 
spontaneous emission, the QS method seizes to work, as the /3^ = condition fixes a k 
to be disproportional to identity 



39 . 



4-1.3. Channel Extension method 

The channel extension (CE) method of 39] not only applies to the widest class of 



quantum maps containing all <£>-non-extremal ones, but also provides more stringent 
bounds than the CS, RLD and QS equivalents, as respectively proven in (39[, 
Appendix C|and above. In this method, see Figure [TJd, each single channel is extended 



by an auxiliary ancilla as introduced in Section l3~3l In |35|, it has been proven that one 
can then effectively bound the A^-channel QFI (]T5]) via the N- extended- channel QFI, so 
that 

j[AH<j[(A^ir] < 4 min {TV ||^|| + iV(iV - 1) ||/3^|| 2 } , (23) 
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Table 2. Quantum phase estimation precision enhancement from the CE 
method. For all noise models specified in |Appendix A| the asymptotic CE bounds 
on the maximal quantum precision enhancement factors, %[•] = •y/.FasM /~F[*], are 
presented. For a general quantum map A v , the CE bound only upperdimits x[A v ] 
as ^asfAy,] < J 7 ^- Yet, for dephasing and lossy interferometry, as indicated by "=", 
the corresponding values of x[A v ] have been shown to be attainable 7^, ZSl- For an 
extended channel, xl-^-v ® 2] is determined by the CE bound as 'F as [A l p (g> I] = J^ E . 



where again ctK = Y^i=i K K i{}P) ■> ^k — ^Y11=i Ki(<fi) an d h i s the generator of 

local Kraus representation rotations Crucially, if there exists a Kraus representation 
for which the second term in (1231) vanishes, ^[A®^] must asymptotically scale linearly 
in N. This requirement corresponds to the constraint Pk = already accounted in the 
QS method, which for any linearly independent Kraus operators {Ki} r i=1 is equivalent 
to the existence of h such that l35| 



2J = 1 



(24) 



What is more, for any channel that admits an h fulfilling (]24j) . one can show basing on 
the results of [35j that the second inequality in (1231 is saturated in the N— >oo limit, so 
that the asymptotic extended channel QFI then reads 



T 



lim - 

N-too 



(A v ® X) 



iV 



4 min 

h 



(25) 



Importantly, (1251) becomes the required asymptotic bound X"k° und °f which we 

refer to as the asymptotic CE bound — X\^ E . We explicitly wrote the full form of 
( [25]) to emphasize its similarity to the extended single channel QFI (Tl2|) . The essential 
difference in (125]) is the constraint (T24l) yielding F^K^ ® X] ^X'fA^ (g) X] and dictating 
the maximal quantum precision enhancement for an extended channel: 



XI 



lim 

iV->-oo 



~ext 
cl 



T\lA v ®X\ 



> 1. 



(26) 



Similarly to (1T21) . (125]) is computable by means of semi-definite programming |39|, so 
that one can efficiently determine both (1181) and (T2"S"]) . The CE-based bounds on xfA^] 
and the exact values of x[A^ <8> X] for the relevant noise models are presented in Table [2] 
Although generally the CE method only upper-limits the maximal quantum precision 
enhancement (TTT]) . it has been proven to quantify x[A v ] exactly in the case of dephasing 
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(72| and lossy interferometer channels (73j|. This has been achieved by showing that input 
states utilizing spin- and light- squeezing respectively yield a quantum enhancement that 
asymptotically attains the corresponding CE-based bounds presented in Table [2j On 
the other hand, as indicated in Table [U these channels are also examples of ones for 
which the extension does not improve the precision at the single channel level, so that 
X[A<p] =x[^tp ® m Table [2j The question — when the lack of precision improvement 
due to extension at the single channel level translates to the asymptotic regime, i.e 
J-'fAJ =J r [A ip <g) X] -<=>- X'asfA^] =J r as [A !f (g) X], we leave open for future research. 



4-2. Finite-N Channel Extension method 



In Section 14. 1[ we have presented the CE method as the most effective one out of all 
discussed that provides the tightest upper limits on the maximal possible asymptotic 
quantum precision enhancement. On the other hand, in the case of experiments such 
as optical interferometry with moderate numbers of photons 0, 0|, the asymptotic CE 
bounds, despite still being valid, are far too weak to be useful. For very low values of 
N, the precision can be quantified numerically, for instance by brute-force type methods 
computing explicitly the QFI. However, in the intermediate N regime — being beyond 
the reach of computational power, yet with iV too low for the asymptotic methods to 
be effective — more accurate bounds should play an important role. 

We propose the finite-N channel extension method, which despite being based on 
the properties of a single channel, still provides bounds on precision that are relevant in 
the intermediate N regime. We utilize the upper-limit (123|) on the TV-extended-channel 
QFI and construct the finite-N CE bound, J 7 ^, that reads 

< 7^ E = 4 minjlKH + (N - 1) \\f3 R \\ 2 }. (27) 



N ~ h 



Following the suggestion of 39j, in contrast to the asymptotic CE bound defined 
in (j2"5|) . we do not impose the SQL-bounding condition 0^ — f[2"4")) . but we seek at each 
N for the minimal Kraus representation that is generated by some optimal h = h opt (N) 



being now not only channel but also iV-dependent. Still, as shown in Appendix E[ T 1 



rCE 

N 

can always be efficiently evaluated numerically by recasting the minimization over h in 
(|27j) into a semi-definite programming task. Moreover, as the finite- N CE bound varies 
smoothly between N = 1 and N = 00, at which it is tight, it provides more accurate 
bounds than its asymptotic version. For N = 1, coincides with the extended channel 
QFI flTJ) — J = n=i=HK® 1 ^ whereas in the N — » 00 limit it attains the asymptotic 
CE bound m~^«> = ^- 

What is more, when considering channels for which the asymptotic CE method 
fails, as it is not possible to set (3^ = in (123|) for any Kraus representation, (l27j) still 
applies and it is the finite- N CE method that provides the correct CE-based bound in 
the A^— j-oo limit that in principle may then surpass the SQL-like scaling. On the other 
hand, when dealing with estimation schemes in which one can moderate the effective 
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Figure 3. Phase estimation CE method based bounds on precision 

(a) Dephasing: Finite-TV (dashed) and the asymptotic CE bounds (dotted) on 
estimation uncertainty as compared with the precision achieved by utilizing spin- 
squeezed states in a Ramsey spectroscopy setup with probe consisting of N atoms 
experiencing uncorrelated dephasing with 77 = 0.9 (solid). 

(b) Loss: Lossy interferometry with particle survival probability r\ = 0.9, e.g. Mach- 
Zehnder interferometer experiencing photonic loss in both of its arms, with effective 
power transmission 77. The smallest uncertainty in a phase estimation scheme is 
quantified by calculating the QFI for numerically optimized TV-particle input states. 
Again, finite- N (dashed) and asymptotic CE bounds (dotted) on precision are shown 
for comparison. 



amount of loss (i.e. the form of A^,) depending on the number of particles, the asymptotic 
bound may not actually be the tightest within the CE method. The (3%- of 
and (127)1 becomes then a function of N and it may not be asymptotically optimal to 
set it equal to zero by imposing condition (j24p . Yet, the finite- A" CE method being not 
constrained with j3^ = still yields the correct CE-based bound on precision as N — too. 



This fact has been utilized in 40|, where, owing to the A^-dependence of the finite- AT 
CE method provided an asymptotic bound indeed tighter than the naively calculated 
J-^ E . What is more, the F^^oo nas been numerically shown there to be saturable, what 
proves the power of the more agile finite-A" CE method. For phase estimation with 
various decoherence models including dephasing, depolarization, loss and spontaneous 



emission described in detail in |Appendix A[ we observe that the finite- A" CE bound is 



simply related to its asymptotic form as 

where one should substitute for the corresponding asymptotic CE bounds presented 
in Table Q]|J 

For dephasing and lossy interferometry models, we show explicitly in Figure [3] both 
the asymptotic and the finite-A" bounds accompanied by the plots of actual precision 
achievable with explicit estimation strategies optimal in the large A" regime. In the first 



case, we consider a Ramsey spectroscopy setup of 13|, [L4j in which the probe consists 



(t In the case of spontanous emission noise the formula is valid only for N > 2, what we suspect to be 
a consequence of the spontaneous emission channel being an extremal map [60] . 
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of atoms prepared in a spin-squeezed state |72|. The parameter is then encoded in the 
phase of a unitary rotation generated by the total angular momentum of the atoms 
that simultaneously experience uncorrelated dephasing. After measuring probe's total 
angular momentum perpendicular to the one generating the estimated phase change, 
the parameter is reconstructed with uncertainty plotted in Figure [3^. In the second 
lossy interferometry case, each particle is subject to an independent loss process with 
survival probability 77, as in e.g. Mach-Zehnder interferometer with power transmission 



r\ in both arms |29l . l30| . Here the solid line represents the maximal QFI achieved with 
numerically optimised iV-photon states. Although in this case the maximal quantum 
enhancement of Table |2] can also be achieved via an estimation strategy that employs 



squeezed-light as the input with mean number of photons constrained to N = N 73 
we cannot compare its precision with the one bounded through As A r - J-^ E is a 

convex quantity in N, one cannot use it naively to constrain precision after replacing N 
by the mean number of photons N. This contrasts the case of the (constant) asymptotic 
CE bound, which yields N-J 7 ^ being linear in N, so that it also applies to estimation 
strategies employing states of indefinite photon number, as pointed in 39| . 

5. Frequency estimation in atomic models 

We apply the methods discussed above to the case of frequency estimation problems in 
atomic spectroscopy. The general Ramsey spectroscopy setup considered in (I3I Q, 0- 



26l . |40| corresponds to N identical two-level atoms — spin-1/2 systems — with their 
states separated, where typically the detuning u of an external oscillator frequency from 
the atoms transition frequency is to be estimated. We assume that the full experiment 
takes an overall time T, during which the estimation procedure is repeated k = T/t times, 
where t is the evolution duration of each experimental shot. The Quantum Cramer- Rao 
bound (j2J) on precision of the estimate Co can be then conveniently rewritten as 

ACu > . , 29 

where ft[g(t)} = Fq [£>(£)] ft is now the effective QFI per shot duration and (t) denotes 
the final state of the whole probe containing N particles being measured in each shot. 
The total time T plays then the role of k in (j2J) and, after fixing t, the bound (1291) can 
always be saturated as T — > 00. The evolution of the probe can be modelled by the 
master equation of the Lindblad (54| form 



^=Ei|K ) .^(*)]+/: (B) wr(*)], (30) 

n=l 

where 03 is the Pauli operator generating a unitary rotation of the n'th atom around 
the z axis in its Bloch ball representation. The uncorrelated noise is represented by 
the Liouvillian part acting independently on each, here the n'th, particle, so 

that effectively p^(t) = A®£ [ V'La)] with channel A w;t representing the overall single 
particle evolution over time t. To model the decoherence we consider the dephasing, 
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depolarization, loss and spontaneous emission maps, which corresponding Liouvillians 
can be found in Appendix A As the estimated parameter corresponds now to oj = ip/t, 



where <p is the phase of the unitary rotation, the QFI via a parameter change just 
rescales, so that ft[gj\ =Fq[q u1 ] /t = Fq[Q ip ] t. Moreover, due to the commutativity of the 
unitary and the considered decoherence maps, we can without loss of generality utilize 
the results presented for them in the previous sections. Defining the channel QFI for 
frequency estimation tasks similarly to (j7|) as 



f[AJ = max max f t [A u . t [\i/; in )]] 

0<t<T ip- in 



(31) 



we can compute all the appropriate expressions for QFIs and the asymptotic bounds 
of Table [TJ as well as the finite- N bounds of (|28p by substituting for the effective time 
dependence of the decoherence strength r](t), which is determined by the master equation 
( 13U]I (see Appendix A ). Then, any quantity J 7 listed in Table [T] transforms to its t- 
optimized equivalent as f = max Tt. In Table [3] we present the channel QFIs relevant for 

0<t<T 

frequency estimation, their asymptotic and finite- N CE bounds, as well as the maximal 
quantum precision enhancements for each model considered. In the case of dephasing, 
we recover the results of 2J, |38|, |40|, whereas for depolarizing, loss and spontaneous 



emission maps we obtain the QFIs and their bounds, which to our knowledge have 
not been reported in the literature before. However, similarly to the case of quantum 
phase estimation summarized in Table [2j the obtained quantum enhancement factors 
for depolarization and spontaneous emission channels serve only as bounds, as they are 
not guaranteed to be saturable. 



6. Estimation of decoherence strength 

Lastly, we would like to emphasize that the CS, QS and CE methods described in 
Section H] also apply to estimation tasks in which the estimated parameter is not encoded 
in the unitary, noiseless part of the system evolution. An example of such schemes are 
the experimentally motivated ones, in which one tries to quantify the effective strength 
of noise or loss present in the apparatus. That is why, we consider again the channels 



described in Appendix A , but this time with the parameter to be estimated being 
the decoherence strength rj. This kind of problems has been widely considered not 
only in the estimation theory 
discrimination 



4J, l74H76|. but also when examining issues of channel 



45l . |46| with particular application in quantum reading 50M52| . As 



compared to unitary rotations, the nature of the estimated parameter is dramatically 
different. In unitary parameter case, the use of entangled input state of N particles 
results in an effective A^-times higher "angular speed" of rotation leading to the HL in the 
absence of noise. In decoherence strength estimation tasks, a change in the parameter 
value can be geometrically interpreted in the space of all valid quantum channels as 
a "movement" in the direction away from the boundary of the space of the relevant 
CPTP maps, which "speed" cannot be naively amplified A^-times when employing N 
parallel channels. Hence, as in the case of lossy unitary rotation estimation, the optimal 
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Table 3. QFIs, CE bounds and quantum enhancements in frequency 
estimation. In frequency estimation tasks the precision is maximized by adjusting the 
single experimental shot duration t. The i-optimised (extended) channel QFIs as well 
as their finite- N and asymptotic CE bounds are presented, where w x [N] = 1+H^[ :r c ]^ r ] , 

w = 1 + 2 W ' 



2^0 



and W[x] is the Lambert W function. As in the case of 
depolarizing channel not all solutions possess an analytical form, only their numerical 
approximations are shown with a w 2.20 and j3 sa 1.32. Right of the double-line the 
maximal quantum precision enhancements are listed for the maps considered. In the 
case of dephasing noise the ultimate ^fe factor has already been reported in 0,|3^|. For 
unextended depolarization and spontaneous emission maps, the derived enhancement 
factors may possibly be not achievable. 



entangled inputs must lead not to scaling but constant factor quantum enhancements, 
which again can be quantified by the methods of Section HI This also explains that for all 
the four noise models considered, the purely geometrical notion of classical simulability is 
enough to bound most tightly the maximal asymptotic precision of estimation. However, 
as for them also J-^ s = J-^ s = = J 7 [A,, <g)Z], the CS-based asymptotic quantum 
enhancement corresponds to the classical estimation strategy that employs independent 
but extended channels. The fact that factorizable inputs — uncorrelated in between the 
extended N channels but possibly requiring entanglement between each single particle 
and its ancilla — are optimal for noise estimation with extended channels, has already 
been noticed for the low-noise \t\ and generalized Pauli |78| channels, of which the 
latter contain the dephasing and depolarization maps here studied. 

In the case of dephasing channel, we further realize that the extension at the 
single channel level is also unnecessary, as J 7 [AJ = J r [A J? ®X] = 1/(1 — rj 2 ), and 
the geometrically dictated bound of CS is attainable classically just by employing 
unentangled qubits in any pure state lying on the Bloch sphere equatorial plane. 
Similarly, the spontaneous emission and loss maps also turn out to be fully classical. In 
the first case, the asymptotic CS bound coincides with the extended and unextended 
channel QFIs derived in 79|, whereas for the loss channel we obtain J-^ s = J 7 [AJ = 
T [A,, <S> I] = 1/ (77(1 — 77)), which at the single channel level is achieved by a photon in 
any mixed state with the ancilla being redundant. On the one hand, this emphasizes 
that entanglement between the photons entering the interferometer is unnecessary and 
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Channel considered 






jrCS jrQS jrCE 

flS cLS 3iS 


Dephasing 


i 

i-n 2 


1 

I-77 2 


l^S 


Depolarization 




(1-»j)(1+3tj) Nl 


(1-77x1+377) Elzl 


Loss 


1 

77(1-77) 


1 

77(1-77) 


1 

77(1-77) 


Spontaneous emission 


,(!-,) S 


,(1-,) H 


1 

77(1-77) 



Table 4. Decoherence strength estimation quantified via channel QFIs and 
their asymptotic bounds. Definitions of channels listed in the first column can be 
found in |Appendix A| In contrast to phase estimation examples given in Table [TJ the 
variable to be estimated here is the decoherence parameter 77, (0 < 77 < 1). Due to 
the different nature of the estimated parameter, the geometrical classical simulation 
(CS) method provides bounds that not just most tightly limit the asymptotic extended 
channel QFIs, but actually coincide with its value. The results prove that only in the 
case of depolarization channel the precision can be enhanced with the use of quantum 
estimation strategies, as for all other cases J 7 [A V ]=J 7 [A J] ® I] = J 7 ^. 



agrees 



with the results of , [75| confirming that the total photon number fluctuati 



10ns 



are really the ones that limit the precision. These can be reduced by employing Gaussian 



rn 



states [74J or in principle fully eliminated by the use of Fock states 
CS-bound. 



751 that attain the 



The case of depolarization map is different, as it is known that for qubits |69l . 180 



the precision of estimation may be improved by extending the channel, i.e. J 7 [A v ] < 
T [A r; ® I]. This leaves the space for possible enhancement thanks to the use of 
entangled probes between unextended channels and indeed this fact has been observed 
already when considering two depolarization channels used in parallel 69|. The results 
are summarized in Table |H 



7. Summary 



We would also like to point that the SQL-like bounds, universally valid in practical 
metrological scenarios, allow to avoid some of the controversies characteristic for 
idealized decoherence-free scenarios. When decoherence is not present and the probe 
states with indefinite number of particles are considered, such as e.g. squeezed states of 
light , the exact form of HL needs to be reconsidered 8lN83| since the direct replacement 
of N with mean number of particles N may make the HL invalid. Moreover, the final 
claims on the achievable precision scaling may strongly depend on the form of a priori 



84 



These 



parameter knowledge assumed, and lead to some apparent contradictions 
difficulties do not arise in the realistic metrological schemes, as the asymptotic SQL- 
like bounds are valid also when N is replaced by N for indefinite particle number 
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state 



39 



1861] . The bounds derived in the so-called local approach (small parameter 
fluctuation) based on the calculation of the Quantum Fisher Information (QFI) are 
saturable in a single-shot scenario unlike the decoherence-free case when only after 
some number of independently repeated experiments one may expect to approach the 
theoretical limits j8?]|. This is due to the fact that by employing input states of grouped 
particles, which possess no correlations in between the groupings, and by letting the 
groups to be of finite but sufficiently large size, one can attain the ultimate asymptotic 
SQL bound up to any precision. Since saturability of the quantum Fisher information 
bounds for independently prepared probes is well established 6lM63ll88|. the operational 
meaning of the QFI is clear also in the single shot scenario. 

We also conjecture from the point of view of the asymptotic SQL-like bounds that 
the specific form of the a priori knowledge should not play an important role. In 
particular, we expect that various methods such as Bayesian 89M92| or information 
theoretic 93|, |94| should recover the bounds compatible with the ones obtained via the 
local approach considered in this paper. This statement is known to hold in the case 
of optical interferometry with losses 3l| , |32| and it is an intriguing question whether 
analogous claims can be made in more general scenarios. 
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Appendix A. Channels considered 

We adopt the standard notation in which 1^ represents a.dxd identity matrix and {<Jj}^ =1 
are the Pauli operators. In Section fl~2l parameter ip to be estimated is the rotation angle 
around the z axis of the Bloch ball generated by the unitary operator U ifi = exp[io~3<p/2]. 
We consider maps T> n with decoherence parameter r\ that commute with such rotation, 
whence A v [g] = V^^U^gU^ = U v V v [g]U^, and are defined accordingly by the Kraus 
operators presented below. For each case, we also specify the purification determining 
the extended channel QFI ffT2l) (X'fA^ £§) X] in Table [T|) in the form of the optimal 
generator of Kraus representation rotation h, as introduced in fllOp . Dealing with 
frequency estimation tasks discussed in Section [5] we construct the effective one-particle 
Kraus operators by substituting (p—>ut and r)—>ri(t) in the nominal ones, where u is the 
estimated frequency detuning. For all models, we explicitly write the Liouvillian 
determining the noise affecting each particle, see fl3U]) . and the effective form of rj(t). 
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When discussing decoherence strength estimation in Section [6J we consider solely each 
of the following noise maps with rj being now the parameter to be estimated: A ip=v = T> ri . 



Appendix A. 1. Dephasing 

• Decoherence parameter, 77, represents the final equatorial radius of the Bloch ball 
shrank uniformly in the xy plane by the channel. 

• Kraus operators: 



I + 77 



I 2 , K 9 



1 — T) 



O-3 ■ 



2 V 2 

Optimal purification determining the extended channel QFI (fT2l : 



it 



(A.l) 

(A.2) 
(A.3) 



h=— r -a 1 . 

• One-particle Liouvillian for frequency estimation tasks: 

Appendix A.2. Depolarization 

• Decoherence parameter, rj, represents the final radius of the Bloch ball shrank 
isotropically by the channel. 

• Kraus operators: 



I + 377 



2 • 



Ki 



1 — 77 



O-i-l 



(A.4) 



t=2...4 



Optimal purification determining the extended channel QFI (TT2T) : 

/o 00 A 



h 








with ^ 



V / (l + 3r / )(l-r / ) 



1+7/ 



(A.5) 



\f 0/ 

One-particle Liouvillian for frequency estimation tasks: 



.-. 7/(t) = e 



(A.6) 



i=i 



Appendix A.3. Loss 



Decoherence parameter, 77, represents survival probability of each of the particles 
that are subject to independent loss processes. The channel on a single particle 
is formally a map from a two- to a three-dimensional system with output's third 
dimension corresponding to the vacuum mode responsible for the particle loss. 
Although, in this case one should strictly write K v , = T> rj \U ip QU^\ = U (p T> r) [g] with 
Up acting on a different Hilbert space, the effects of U v and U v are physically 
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indistinguishable, as the particle losses commute with the acquired phase (for 



instance see 



30||). 



Kraus operators: 

/0 \ 











/ o o\ 



\0 VT^jJ 

5t 
/ 

V 



\V^v o J 



y/Tj 





(A.7) 



h 



0"3 



Optimal purification determining the extended channel QFI (TT2l ) : 

\ 


j 

One-particle Liouvillian for the frequency estimation tasks: 



(A.8) 



m=0 ^ 



(n) (n) £ / (n) (n) 



-7* 



(A.9) 



where cr^ + = |vac)(m| are the generators of transition to the vacuum mode from 
qubit basis states |0) and |1), such that cr ( ^_ = a ( ^ > ]_. 

Appendix A. 4- Spontaneous emission (amplitude damping) 

• Decoherence parameter, rj, represents the radius of the disk obtained by projecting 
the deformed Bloch ball outputted by the channel onto the xy plane. 

• Kraus operators: 

11^ „ (0 VT^v' 



(A.IO) 



y/fj I 

Optimal purification determining the extended channel QFI (TT2l ) : 

k= l(« °) with f =i^. 

One-particle Liouvillian for the frequency estimation tasks (a± = | (<Ji ± i^)): 

M-±{a?V?, e }\ .: v{t) = <r*. (A.12) 



(A.H) 



Appendix B. Equivalence of RLD-based bound applicability and local 
classical simulability of a channel 

Given a channel — a CPTP map A^ : H- m — > "H out — we define its Choi-Jamiolkowski 
(C-J) matrix representation 60j as Q\ = A^ ® X[|I)(I|] = Y2i \Ki{^)){Ki{ip)\, where 
{Ki((f)} r i=1 are r linearly independent Kraus operators of A^, and adopt a concise 
notation for bipartite states, in which \<p) = XVj^i^" ^\4>\j) V) |j) = ® I |I) =1 ® T |I) 
with 1 1) = Yli=x' Hm \^) I*)- For simplicity, from now onwards we drop the explicit cp 
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dependence of operators, assuming that the estimation is performed for small variations 
5ip around a given, fixed <p. In Sup. Mat. of 39| (Equation (S9)) it has been proven 
that the condition for any channel to be (^-non-extremal (classically simulable) at ip is 
equivalent to the statement that there exists a non-zero Hermitian matrix such that 



On the other hand, the RLD-based bound exists there if and only if 37 



Pa , Cl a Pa 







(B.2) 



where Pq ± is the projection onto the null-space Q±, i.e. the subspace orthogonal to Q/\ , 
so that Vj : P^ ± \Ki)=0. The (1B.1|) implies (IB. 21) . as by substitution 

P n± ($>v IKMKjh Pn ± = £ fc> p (K P \K P ) fi p )j P n± |^)(^-| P n± = , (B.3) 

\ ij / ij \ V / 

thus any y?-non-extremal channel must admit an RLD-based bound on its extended QFI. 
In order to prove the other direction, we split the derivatives of each C-J eigenvector 



into components supported by tt\ and in the null-space Q±: 
Hence, after substituting for Cl^ the (IB.2[) then simplifies to 



^>=E i ^l«i>+|^ L >- 



(B.4) 



and since A^A = implies A* = A = and {\K i )} i are orthogonal, we conclude that all 
= 0. Thus, (IB. 21) implies that Ki \ = X] 7 - v ij which due to the local ambiguity 



of Kraus representations ([9]) is equivalent to 



j 



its Hermitian and anti-Hermitian parts z/ = z/ H + z/ AH , so that 
^ H ^0 for any non-trivial channel. Finally, we can write 



h. Therefore, without loss of generality, we may set h = — |^ AH after splitting v into 



^Avij — i hij) Kj \ for any Hermitian 



y .z/ H 



Ki > with 



(B.5) 



ij 



and satisfy the condition (IB.lj) .! 



Appendix C. RLD-based bound as a special case of asymptotic CE bound 



For a channel that admits an RLD-based bound, in order to obtain the CS condition 
(IB.5P in Appendix B, we chose h = —^i ,AR that actually satisfies the /3^ = constraint 



(I24p of the CE method. This can be verified by taking the Tr% out {. . .} of the both sides 
of the identity 

\K 3 )(IQ\ = £i (u v - 4) \K s ){Ki\ = l -Y,\ k ){Ki\ - |/^)(^|, (Cl) 



ij 



•j 
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which results in (I24p . This is consistent, as the CE method must apply to any (^-non- 
extremal channel |39j — admitting an RLD-based bound. Furthermore, the asymptotic 
CE bound (125]) is at least as tight as the RLD-based bound (TT3j) on the extended channel 
QFI (P2D. We prove this by substituting (IB. 51) into the definition of J 7 ^ [A^ <g> X] in 
(CGI, so that 



J- RLD [A V 



ij 



(C.2) 



n- 1 \K P 



5 



Hence, J rRXD [A ¥ , (gi 2] is 
an example of an asymptotic CE-based bound with a possibly sub-optimal Kraus 



where we have used the fact that \^Kj 



representation chosen such that Vj 



K, 



Kj ) and (3 



A 



0.1 



Appendix D. Optimal local QS of a channel 

A channel A^ of rank r, in order to be locally quantum simulable within small deviations 
Sip from a given tp, must fulfil the condition (see Section 14. 1 . 2j) 

r'>r 

AM = Tr E$ E CT {t/ (q ® I^X^I) tf} + 0(5<f 2 ) = KM gKi&y + 0{5ip 2 ), (D.l) 

i=l 

where Ki(tp) = (i\ U \ifjp) and {|i)}[ =1 form any basis in the r' dimensional Ke 4> x "He ct 
space containing ip^. Hence, A v must admit a Kraus representation {Ki}l' =1 (with 
possibly linearly dependent Kraus operators, as for generality we assume r' > r) that 
coincides with the one of (ID. ID up to 0(d~ip 2 ), i.e. satisfies K^ — Ki and Ki = Ki for all 
i. We construct a valid decomposition of |^} into its (normalized) components parallel 
and perpendicular to ip, 



^> v y=i a|Vv) — 1 H^)' wriere we can choose a, 6g1R because 
of ^(Vvl^) = and the irrelevance of the global phase. Then, the asymptotic bound 
jrbound of (Tjgj) determined by the local QS (lD~Tj) at y? simply reads F Q [|^)] =46 2 and 
the required Kraus operators {Ki}l' =1 of A^ must fulfil conditions Ki = (i\U \ip v ) and 
Ki = (i\ U ip (p )=iaK i —ib (i\ U |"0i")- Hence, for the local QS of channel A v to be valid, 



we must be able always to construct 



U 



k 3 f^ + i^ 



(D.2) 



with first two columns fixed to give for all i the correct (i\U and (i\ U |^) 
respectively. Due to locality, all entries marked with • in (1D.2|) can be chosen freely 
to satisfy the unitarity condition U'U — UU^ =1. Yet, this constraint still requires the 

Kraus operators to simultaneously fulfil i Y^i=\^ i^-i = al and X^=i-^i-^ = + fl2 ) ^- 
Without loss of generality, we may shift their phase at ip, so that Ki — > e~ mip Ki and 
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the conditions become independent of a, i.e. ^YH=i^i^i = an d YH=i^i^-i = 
Furthermore, these constraints do not require r'>r, as rewriting for example the first 
one as i^i=i\VV U N) (^1 ^MVv) = 0' one can always resolve the identity with some 

basis vectors 5^i=il0 (*l = Si=il e «) ( e «l an d define linear independent Kraus operators 
{Ki = (ej| U \ipip)Y- =1 also fulfilling the necessary requirements. 

Finally, we may conclude that is locally quantum simulable at <p, if it admits 
there a Kraus representation satisfying conditions (|2"T|) stated in the main text, which 
due to locality can be generated via by some Hermitian rxr matrix h.M 



Appendix E. Finite- N CE method as a semi-definite programming task 



The finite-N CE bound has been defined in (1271) as 



CE 
N 



4 min 1 1 



a 



K 



(N 



K 



(E.l) 



where ||-|| denotes the operator norm, = ^2iK { Ki and (3^- = i^ i KiKi- Given a 
channel A^ from a <ii n - to a <i out -dimensional Hilbert space and the set of its linearly 
independent Kraus operators (fi utX^m matrices) {Ki} r i=1 , in order to compute J 7 ^, we 
should minimize (IE. II) over locally equivalent Kraus representations of A^ generated 
by all Hermitian, rxr matrices h. 



Basing on the results of |39j, where the = constraint (124)) is also imposed on 
(lETTj) , we show that can always be evaluated by means of semi-definite programming 
(SDP). Adopting a concise notation in which K is a column vector containing the 
starting Kraus operators Ki as its elements, we can associate all locally equivalent Kraus 
representations K in (IE. ID with those generated by any h via K = K and K = K — i/iK. 
By constructing matrices (Id represents a dxd identity matrix) 



A 



K 



~ t 
K 



B 



iKK 



iK f K 



(E.2) 



which positive semi-definiteness conditions correspond respectively to 



Oi i 



KK < Xnh 



we rewrite (IE. II) into the required SDP form 
J-£ E = 4min{A a + (iV-l)A 6 }, 



K f KKK < A b Id in , 



(E.3) 



(E.4) 



subject to: A > 0, B > 0. 



For the purpose of this paper we have implemented a semi-definite program using 
the CVX package for Matlab 95|, which efficiently evaluates ( 1E.4|) given the set of Kraus 
operators and their derivatives of a generic channel A v . The fact that only K and K 
are involved in (IE.4[) is a consequence of the QFI, and hence the bound J 7 ^, being a 
local quantity. 
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Lastly, one should note that by slightly modifying the program in (IE. 4}) we are able 
to also efficiently evaluate: the extended channel QFI (lT2"j) . as J-^A^ <8>I] = J 7 ^; and 
the asymptotic extended channel QFI (125|) . J^fAy, ®X] = , by setting = 1 and 



imposing the /3r- = constraint (T2"4")) as already pursued in 39l 
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